Conductance of a quantum dot in the Kondo regime connected to dirty wires 
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We study the transport behavior induced by a small bias voltage through a quantum dot connected 
to one-channel disordered wires by means of a quantum Monte Carlo method. We model the 
quantum dot by the Hubbard- Anderson impurity and the wires by the one-dimensional Anderson 
model with diagonal disorder within a length. We present a complete description of the probability 
distribution function of the conductance within the Kondo regime. 
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The Kondo effect is one of the most paradigmatic phe- 
nomena of electronic correlations. It was proposed long 
time ago to explain the peculiar behavior of the resistiv- 
ity of magnetic impurities in metals as a function of the 
temperature [l[ . Towards the end of the last century, the 
interest on this effect was paramount to the mesoscopic 
community, after it was observed in transport experi- 
ments in quantum dots 0, HJ . 

A quantum dot is a confined structure where electrons 
experience a strong Coulomb repulsion, which is in con- 
tact to metallic wires. For temperatures lower than the 
Kondo temperature Tk, the Coulomb interaction origi- 
nates an effective coupling between the spin of a localized 
electron at the dot and the spin of the electrons of the 
wires. The result is the formation of a resonant singlet 
state which manifests itself as the opening of a transport 
channel, corresponding to a conductance G — e 2 /h for 
each spin component. The electrons of the wires that 
intervene in the formation of these singlets define the 
so called "screening cloud," which extends over a length 
that scales as 1/Tk- The build up of this cloud and 
the development of the Kondo effect are sensitive to fea- 
tures related to the finite size of the wires. The peculiar 
behavior of G for a dot in the Kondo regime when the 
screening cloud does not fit into the finite wires has been 
analyzed in Ref. 0, and the features there predicted were 
later confirmed by a numerical study jB|. These studies 
assume perfect clean wires where the finite-size effects 
are solely defined by constrictions within the wires at a 
finite distance from the quantum dot, but do not analyze 
the effect of disorder, which is an unavoidable ingredient 
in most of the realistic experimental settings. 

In low dimensional conductors, backscattering induced 
by disorder produces localization within a length £ m /j,. 
This is one of the most dramatic consequences of the 
quantum coherence characterizing the transport in meso- 
scopic devices, which enables the interference of the elec- 
tronic wave function. This effect has been the subject of 
many investigations in the framework of non-interacting 
systems [0, llij LLJI > while there are also some studies for 
interacting ones [7]. In the first case, the probability dis- 
tribution function (PDF) for the conductance at temper- 



ature T = 0, P (G°(T = 0)) , is analytically known for ID 
systems where the disorder is described by the disordered 
Anderson model @, @] ■ That description has been also 
extended to the case of finite voltage and finite tempera- 
tures The function P (G°(T)) is completely char- 
acterized by a single parameter (/£, m fp = — (l n G°(0)) ex- 
cept for anomalies at the band edges and the band center 
WL 3], being £ m /j, the mean free path and £ the length 
of the dirty part of the wires. In the context of the Kondo 
effect, disorder in the environment of the interacting im- 
purity is expected to affect the development of the Kondo 
resonance. This problem has been analyzed in the liter- 
ature under the name of the "Kondo box," where the 
probability distribution for the Kondo temperature Tk 
has been shown to depend on the dimensionality and the 
size of the box [H-[l8l. 

The interplay between the electronic correlations tak- 
ing place in the Kondo regime, and the localization in- 
duced by the disorder has not been so far considered in 
the context of the transport properties. In the present 
work, we precisely address this important aspect. We 
identify a crossover in the behavior of P (G(T)) as the 
length of the dirty piece of the wire increases over the 
mean free path £ m / p and as the temperature overcomes 

We describe the quantum dot by a Hubbard- Anderson 
model connected to left (L) and right (R) wires, which 
are dirty within a finite length. The ensuing Hamiltonian 
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where the Hamiltonian for the dot includes the effect of 
the Coulomb repulsion U and the voltage gate V g . It 
reads 



H d =V a 



J2 n ^ 

=r=t4 



Un d ,fnd,. 
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Each wire is modeled as a one-dimensional Anderson 
Hamiltonian with diagonal disorder within a length 1/2 = 
Na, being N the number of sites of a tight-binding lattice 
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of lattice constant a, which we set to unit. The Hamilto- 
nian has the form H a = H a , w + H a c + H a ^ ri being the 
dirty piece 

N-X N 
H a>w = —t \ c i,g t a c j+\,a,a +H.C.] + ^~~^ e j,a c ^j.a,o c j,a,<7i 
3=1 3=1 

(3) 

where the local energies e^ a are randomly distributed 
within a range [— W, +W] for 1 < j < N. The disordered 
chain is connected through H a c = — t c>a [c' N a CN+x, a + 
H.c] to a clean semi- infinite chain (the reservoir) H a r — 



j=N+l \yj,a,rr 



[Cj a a Cj+i^ a ^ a + H.c.]. The Hamiltonian de- 



scribing the contact between the dot and the wires reads 
H cont = -t i c \, a ,,da + H.c.}. (4) 

a=L,R 

The conductance through the dot for each spin chan- 
nel, in units where e = h = kg = L for an infinitesimal 
applied voltage V, such that the chemical potentials of 
the left (L) and right (R) reservoirs are pl,r = A* ± T^/2 
reads El: 



G 



(5) 



where /(w) = [1 + e' 3 ^ ^] 1 is the Fermi function, 

T(u) = r £ (w)r fl (w)/[ri(u) + r B (w)], with r a (w) = 

— 2t 2 Im[g^(uj)], where (w) is the retarded Green func- 
tion of the wire a connected to the reservoir but uncou- 
pled from the dot. Lastly, p„{u,T) = -2Im[G| CT (w)], 
is the local density of states (LDS) at the quantum dot. 
While the non-interacting Green functions can be rather 
straightforwardly evaluated from a recursive procedure, 
the retarded Green function G§ a (uj) depends on the in- 
teractions and, thus, on the temperature. 

Quantum Monte Carlo methods are powerful tech- 
niques to study the effect of interactions in quantum 
transport 
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For a given disorder realization, we 
evaluate the Matsubara Green function by means of the 
quantum Monte Carlo method of Ref. [2l| and we then 
use a polynomial fit to compute this function for real 
lu. This procedure is very precise within the low-energy 
range, \u\ < U. We evaluate the LDS p a {u,T) and use 
this function to compute G(T). We repeat this procedure 
for up to 5000 realizations of disorder. 

In the limit of vanishing diagonal disorder e;. Q 
0. the model reduces to the Kondo impurity with 
an effective exchange constant J = 2t 2 [l/(U + 
V g ) — 1/Vg] [H, which displays a resonance at the 
Fermi level below the Kondo temperature defined from 
2 = J J*^° dwtaxih.[u)/(2TK)]p w (u))/u), where p w (uj) = 

~ I m E Q 3a( w )] i s tne density of states of the wires. 
This feature has a significant impact on the transport 
properties. In a clean system without mismatching 
(t c .a = t), the solution of this equation is and the 



conductance per spin is exactly the conductance quan- 
tum for T < Tj(. The effect of a mismatching in the 
contact between the wires and the reservoirs, introduces 
finite-size features in p w (uj), thus affecting the behavior 
of the conductance. These effects become relevant within 
a temperature range of the order of the mean level spac- 
ing of the wires when this scale is larger than T^, which 
corresponds to short enough wires (3, IE]. 

For dirty wires, a given disorder realization ei tU ^ 
defines the behavior of p w {u) which in turn determines 
the corresponding value of Tk- For an ensemble of dis- 
order realizations, there is a distribution of Kondo tem- 
peratures, which in ID has a crossover from a log-normal 
distribution centered at T^-, to a distribution with a long 
tail at low temperatures as the disorder amplitude W 
increases, while keeping constant the size of the dirty 
18[. We have verified that the same behavior 
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piece 

is obtained when W is kept fixed and the length of the 
dirty wires increases. In this case, the crossover takes 
place as the length of the wires overcomes the mean free 
path £ m / p . The mean of the distribution, however, re- 
mains approximately T^. In what follows, we analyze 
the consequence of that crossover in the behavior of the 
conductance of this system. For simplicity, we fix p = 0. 
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FIG. 1: (Color online) Top panel: Mean value of the con- 
ductance (G(T)} as a ruction of the temperature T for wires 
of different lengths £ = 14, 40, 100, 200 and 400. The corre- 
sponding values for a non-interacting dot are shown in dashed 
lines for comparison. Also for comparison, we present in full 
line the behavior of the conductance as a function of T for 
a dot connected to perfect clean wires. Bottom panel: Vari- 
ance C2(T) for the same distribution as in the top panel. The 
Coulomb interaction is U = 8t and the disorder strength is 
W — 2t. The mean free path of these wires in the non- 
interacting limit is £ m /p = 144. 

The most interesting features introduced by the Kondo 
effect are expected to emerge at finite temperature. In 
Fig. [1] we show the mean value of the PDF of the con- 
ductance (G(T)) as a function of T for disordered wires 
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of different lengths connected to an interacting dot in 
the Kondo regime, along with the corresponding variance 
cr 2 (T) = ((G(T) - (G(T))) 2 ). The most dramatic feature 
noticed in this figure is the exponential drop of the mean 
conductance (G(T)) as the temperature grows above the 
Kondo temperature T^. This drop is very pronounced 
for short disordered chains with t < £ m / p , and becomes 
milder as the length overcomes £m/ P ; hi which cases the 
conductance is significantly smaller than the conductance 
quantum at T = 0. This is in strike contrast with the 
behavior of the mean conductance for U — 0, (G°(T)), 
also shown in the Fig. Q] (see dashed lines), which re- 
mains approximately constant within a range of temper- 
atures of the order of the bandwidth of the reservoirs 
12! . The high temperature drop of the conductance, 
can be analyzed taking as a reference the behavior of 
the conductance of a dot with the same parameters but 
connected to clean wires with perfect matching with the 
reservoirs, which is also shown in full line in the Fig[TJ In 
the clean system, it is known that the process of singlet 
formation with one spin at the dot and one spin at the 
wires producing the Kondo resonance is effective only at 
low temperatures T < T%. Above the Kondo temper- 
ature the dominating process is the Coulomb blockade, 
which is a consequence of the high energy cost to intro- 
duce an additional electron in the dot, once it has been 
previously occupied by another one. In this regime, the 
Kondo peak at the Fermi energy of the LDS at the dot 
melts into a valley which lies between two peaks sepa- 
rated by U [l|, and the consequence is a dramatic drop 
in the conductance. For dirty wires in which case we 
have a distribution of Kondo temperatures with mean 
Tj^, we can also identify the drop in the conductance as 
a temperature-driven transition from the Kondo to the 
Coulomb blockade regimes. 

A more subtle feature to analyze is the difference be- 
tween the interacting and non-interacting PDF, which 
takes place for short wires within the low temperature 
regime, T < T^. Although small, this difference can be 
appreciated in the Fig. [T] at the level of the mean value 
but also at the level of the variance 02 (T) in relation 
to the corresponding non-interacting value cj^T). Re- 
markably, for short wires, (G(T)) > (G°(T)}, while the 
variance for the interacting case 02 (T) is a decreasing 
function of T, as the non-interacting one, although sys- 
tematically smaller. 

In Fig. [2] we focus on the behavior of the variance 
of both an interacting and a non-interacting dot at a 
fixed temperature T < as a function of the length 
I. It is important to mention that the temperature is 
low enough to allow us to express the conductance by 
G{T) = T(0)p a (0,T), i.e. to disregard the corrections 
of approximating the function df /du> ~ 5(iS). It is clear 
that the low temperature difference between the interact- 
ing and non-interacting variance of the PDF persists up 
to the mean free path £ m / p . Actually, it can be seen that 
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FIG. 2: (Color online) Variance for the interacting (<72(T)) 
and non- interacting (0-2 (T)) system a s a function of the 
length £ for wires at low temperature T « T^-/3 and 2T^/3. 
Other parameters are the same as in Fig. [1] 



for t > £ m / p , the full conductance distribution function 
exactly coincides with the non-interacting log-normal dis- 
tribution, P(G(T)) = P (G°(T)). 
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FIG. 3: (Color online) PDF for the LDS, P(p CT (0)) at T = 
Tx/3 of an interacting quantum dot in the Kondo regime, 
corresponding to three different lengths of the dirty wires i = 
14, 40 and 100. For comparison, we also show in dashdot lines 
the function P(po(0)) corresponding to the non-interacting 
system, as well as in dashed lines the PDF corresponding to 
the Fermi-liquid description. Other parameters are the same 
as in Figs. [T] and [2] 

In order to gain insig ht on the low T < T£ behav- 
ior of P (G(T)) as a function of £, we analyze the be- 
havior of the PDF for the LDS at the Fermi energy 
of the dot, P (p a (0, T)). To this end we resort to the 
Fermi liquid description, expected to be valid within 
this regime. This implies that for low u> and T the 
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self-energy representing the many-body interactions at 
the dot, E a (u,T) = £^(w,T) +iE^(w,T) ) behaves as 
E^(w,T) cx w 2 +tt 2 T 2 [22]. In that limit, the Green 
function for the impurity can be approximated by 



z{T) 



lo - z(T)V g - iz{T) (££(w) + E£(w, T)) ' 



(6) 

where z is the quasiparticle weight, z (T) = 1 — 
9 w E' a (0,T), and = V g - E(,(0) - E^Q.T), being 
E (w) = i 2 E a Sa( w )' ^ n wna t follows, we focus on 
Vg = {7/2 (Vg w 0) which corresponds to the particle- 
hole symmetric point of the model, at which the Kondo 
effect is maximum. From this expression it can be verified 
that, up to 0(T 2 ), the Fermi liquid LDS at the quantum 
dot can be written as 



pf(0,T)=p°(0) 



1 



sg(o) 



(7) 



being pq(uj) the non-interacting LDS, while £^(0,T) rep- 
resents the inelastic scattering effects introduced by the 
Coulomb interaction that take place at finite T. 

The Monte Carlo procedure we follow allows us to 
evaluate, for a given disorder realization, p a (ui,T), while 
p°(0) can be obtained independently by a numerically 
exact recursive procedure. In Fig. 3 we show the his- 
tograms corresponding to these two LDS and it is clear 
that these two distributions arc different for short wires 
with £ < £ m fp, which justifies the difference in the first 
two moments of P(G(T)). The quantum Monte Carlo 
data also provides the distribution for E CT (0,T). The 
imaginary part of this self-energy is peaked around its 
mean value, which we consider along with the data for 
p°(0) and E to evaluate the Fermi liquid density of states 
given by Eq. ([7]). The corresponding histograms, pre- 
sented in Fig. 3, show that P (p^ 1 (0, T)J reproduces 
very well the exact P (p a {Q,T)) within the whole range 
of lengths. The fact that for sufficiently long chains the 
distribution function P (p a (0,T)) along with the Fermi 
liquid approximation exactly coincides with the corre- 
sponding distribution for the non-interacting case can be 
understood within the Fermi-liquid description. Indeed, 
an analysis of the mean value of E^(0,T) as a function 
of the length indicates that this quantity decreases lin- 
early as £ increases and falls to zero at £ ~ £ m / p . This 
is consistent with a picture where the inelastic scattering 
processes accounted by E" are determined by the square 
of the available phase space, which is oc T 2 in a clean 
system and reduces as the system becomes dominated 
by the disorder and localizes. 

Substituting the Fermi liquid density of states (J7J in 
([5]) defines a probability distribution which reproduces 
the corresponding function for the conductance P(G(T)) 
for T < Tk- The particular case of T = is rather 
straightforward after noticing that p& (0, 0) = p°(0), 
which means that the conductance distribution for the 



interacting dot with dirty wires in the Kondo regime 
exactly coincides in this limit with the corresponding 
to the non-interacting case, which is analytically known 

Hi mm. 

To conclude, the finite size of the wires has been iden- 
tified as a source of anomalous behavior in the trans- 
port properties of quantum dots, particularly within the 
Kondo regime 4.5]. In the present contribution we have 
added the ingredient of disorder to that scenario. We 
were able to propose a full description for the PDF of the 
conductance of a quantum dot connected to disordered 
wires within the Kondo regime. For T = this function 
coincides with the one for non-interacting systems. For 
< T < and £ < £ m f p , it can be described by a Fermi 
liquid density of states corresponding to Eq. (J7J , while for 
£ > £mfp there is a crossover to a regime where the lo- 
calization dominates over the inelastic scattering effects 
originated in the interactions and the PDF corresponds 
to the same log-normal function of the non-interacting 
system. For increasing temperatures T > Tj^, the ex- 
ponential drop of the conductance that characterizes the 
Kondo regime is also observed in these systems. The drop 
is milder as the length of the dirty wires increases, since 
the corresponding values for the low-temperature regime 
are significantly smaller than the conductance quantum. 
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and UBACYT, Argentina. 
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